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^ I Abstract. A new way to study the multiplicity production in high energy processes is proposed. It 

Q , is based on the multiplicity distribution description by the schemes taking into account the formation 

■ of quark-gluon system and hadronization. This investigations revealed an active role of gluons and 

' the universal mechanism of their hadronization. 

a^ : 

^ : INTRODUCTION 

T— I ■ 

T— I ' At present there is neither single consistent theory nor convincing model which could 

^ . explain all the results obtained at RHIC and SPS Jil]. It is evident that in these exper- 

ly-^ \ iments the local thermal equilibrium is reached and thermalized matter is produced. 

O ' Thermal statistical models based on this approach explain yields of different hadrons at 

^ ■ high energy [1]. At the same time these models describe well the particle yields in e+e^ 

Ph! and pp interactions where, as generally accepted, thermalization is impossible [2]. So 

p the intensive study and insight of these processes remain to be very importance. 

^ ■ To investigate multiparticle production (MP) at high energy in these processes, a two 

• • ! stage model was proposed 10,0]. It is based on the use of QCD and the phenomenolog- 

. ^ \ ical scheme of hadronization. The model describes well multiplicity distributions (MD) 

^ ■ and their moments in e+e^ annihilation, pp and pp interactions and gives complemen- 

. tary information to a better understanding of MP of relativistic heavy ion collisions. 

■ " " ' This model confirms the fragmentation mechanism of hadronization in e+e annihila- 

tion and the transition to recombination mechanism in hadron and nucleus interactions. 
It explains the shoulder structure in MD at higher energies and the behavior of fi vn pp 
annihilation at few tens GeV/c by including intermediate quark topologies. The mech- 
anism of the soft photons (SP) production as a sign of hadronization and estimates of 
their emission region size is proposed 

The e+e annihilation is the most suitable process to study MP. It can be realized 
through the formation of virtual 7 or Z^-boson which then decays into two quarks: 
e+e" — i> {TP /y) — > qq. The e+e^-reaction is simple for analysis, as the produced state 
is pure qq. To describe the process of parton fission (quarks and gluons) at big virtuality, 
pQCD may be applied. This stage can be named as the stage of cascade. When partons 
get small virtuality, pQCD cannot be applied. Therefore phenomenological models are 
used to describe hadronization in this case. 

The probabilistic nature of parton fission in QCD has been established in |j@]. 



E+E- - ANNIHILATION 



A.Giovannini 0] had proposed to describe quark and gluon jets as markovian branching 
processes with three elementary contributions: gluon fission, quark bremsstrahlung and 
quark pair production. He constructed a system of differential equations for generating 
functions (GF) and obtained solutions of MD for quark jet 
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These MD are known as negative binomial distribution (NBD). The GF for them is 

oo 

Qiq) Y)= z'^'PmiY) = [1 +m/kp{l - z)] ", where m is the mean gluon multiplic- 



ity, 7— QCD evolution variable and kp - NBD parameter. 

In lH MD O was taken for the cascade stage. For the hadronization a sub narrow 
binomial distribution (BD) was added 
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where C^^ - binomial coefficient, n'^ and Np (p = q, g) have a sense of mean multiplicity 
and maximum number of secondary hadrons formed from parton (p) at its passing of the 
second stage. 

We have chosen BD from the analysis of experimental data on e^e^- annihilation. 
Second correlation moments /2 were negative at low energies (less than 9 GeV). The 
choice of such distributions was the only one that could describe the experiment. We 
suppose that the hypothesis of soft decoloration is right. Therefore the hadronization 
stage is added to the parton stage by means of a factorization. We introduce parameter 
a — Ng/Nq to distinguish the hadrons, produced from quark or gluon. MD of hadrons in 
- annihilation can be written as follows (A^^ 
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The expression © describes well the experimental data ^ from 14 to 189 GeV 
ll^ (e.g. Fig. 1). The mean gluon multiplicity m has a tendency to rise, but slower 
than the logarithmic one. Parameter kp was been related with temperature in (Td]. The 
fact that a < 1 is the evidence that hadronization of a gluon is softer than a quark. It 
is surprising that gluon parameters of hadronization {Ng, np remain constant without 

considerable deviations in spite of the indirect finding: A^g ~ 3 — 4 and ~ 1 (e-g- Fig. 
2). Therefore we can draw a conclusion about the universality of gluon hadronization in 
e^e^ annihilation in the sufficient wide energy region. 

It was shown [11] that the ratio of factorial cumulative moments over factorial mo- 
ments changes the sign as a function of the order. The calculation of this ratio by using 
© was done in 10] ■ It has been obtained that the period of oscillations is equal to 2 
before Z*^ -region and increases at higher energies. 




PP INTERACTIONS 

The study of MP in pp interactions is implemented in the framework of the project 
"Thermalization" [12]. This project is aimed at studying the collective behavior of 
secondary particles and advancement to the high multiplicity region (HMR) beyond 
available data id] in pp interactions at 70 GeV/c. The calculation by the MC PHYTHIA 
code has shown that the standard generator predicts a value of the cross section which 
is in a reasonably good agreement with data at small multiplicity n^h < 10, but it 
underestimates the value o{nch) by two orders of the magnitude at Hch = 20. The existing 
models are very much sensitive in this region iHUl . also (Fig. 3). 

We suppose that after the inelastic collision the part of the energy of the initial impact 
protons is transformed to the inside energy. Several quarks and gluons become free 
and form quark-gluon system (QGS). Partons which can produce hadrons are named 
the active ones. Two schemes were proposed |0]. In the first scheme the parton fission 
inside the QGS is taken into account (the scheme with a branch). If we are not interested 
in what is going inside QGS, we come to the scheme without a branch (TSMT). 



At the beginning of research we took a model with active quarks and gluons. Pa- 
rameters of that model had values which differed very much from parameters obtained 
in e+e^- annihilation, especially for hadronization. It was one of the main reasons to 
refuse the scheme with active quarks. So, we reserve quarks remained inside of the lead- 
ing particles. All of the newly bom hadrons were formed by active gluons. That is why 
we began to name this model - the gluon dominance model (GDM). 

The Poisson distribution was chosen as the simplest MD for active gluons which 
appeared for the first time after the collision. The number of these gluons fulfils the role 
of the impact parameter for nucleus. To describe MD in the gluon cascade (fission), Farry 

distribution was used with GF (z) = /m^ [l — z (l — =)\ ^ , where k - the number 
of initial gluons in QGS, m - the mean multiplicity of them in the end of the branch. On 
the second stage some of active gluons can leave QGS ("evaporate") and transform to 
real hadrons. Our BD © for hadrons on the stage of hadronization is as follows 

^>.)^C-(^) (.-^) . (4) 

In © parameter 5 is the ratio of evaporated gluons leaving QGS, to all active gluons. 
From the comparison with data [13] we have obtained that a maximum possible number 
of hadrons from a single gluon looks very much like the number of partons in the glob of 



cold QGP of L.Van Hove lllSll . the branch processes are weak. The fraction of released 
gluons is equal to 5 = 0.47 ±0.01 (the same as in lH^ '). A part of active gluons does not 
convert into hadrons. They stay in QGS and can become sources of soft photons (SP). 

In the scheme without a branch we consider that evaporated gluons have Poisson MD, 
too. Using the idea of the convolution of two stages Pn{s) = Y^PI^{s)PI^ {m) as well as 

m 

the BD for hadronization, we obtain MD in pp-coUisions. The comparison GDM with 
the data [13] (Fig. 3) gives values of parameters Ng = 4.24 ±0. 13, = 1.63 ±0.12 and 
it is in agreement with the values obtained in e+e^ annihilation We are restricted 
in sum to m = 6. Our estimation of the maximal possible observable number of charged 
particles is < 26. 

We have also got MD for neutral mesons and total multiplicity by using mean multi- 



plicities of ;r^-mesons lll7h and active gluons, expected approximate equality of proba- 
bilities of the formation charged and neutral particles from single gluon at hadronization 
and the above-mentioned idea of the convolution |5]. The maximum observable number 
of them is estimated as 16, total - as 42, the parameter = 1.036 ±0.041 (h = TZ^). 

The analysis of the mean multiplicity of Tt*^- mesons versus the number of charged 
particles rich, allows to set limitations to the number of neutral mesons at given rich and 
indicates the absence of AntiCentauro events (Centauro may be in HMR). The obtained 
estimations of probabilities for charged and neutral hadrons production from gluon while 
its passing through hadronization permits to get "the charged hadron/pion" ratio LI 8.1 in 
pp interactions. At 69 GeV/c this ratio is equal to 1.19 ±0.25. 

GDM describes well MD in the region of 100 - 800 GeV/c (Fig. 4 and table 1 in 
101). The number of active gluons, their mean multiplicity, Ng and increase slowly. A 

growth of nll'mpp interactions indicates a possible change mechanism of hadronization 



of gluons in comparison with e+e" annihilation. It is considered that in the last case 
the partons transform to hadrons by the fragmentation mechanism at the absence of 
the thermal medium. Our MD analysis gives n^l ~ 1 for this hadronization The 
recombination is specific for the hadron and nucleus processes. In this situation a lot 
of quark pairs from gluons appear almost simultaneously and recombine to various 
hadrons [19]. The value n'l becomes bigger ~ 2 — 3, that indicates this transition. The 
recombination mechanism provides justification for applying the statistical model to 
describe ratios of hadron yields and the explanation of the collective flow of quarks 
id. 

In O] MP is described by means of a clan mechanism and emphasizes the gluon nature 
of clan. Our GDM allows to give a concrete content for the clan. The clan model uses 
the logarithmic distribution (LD) in a single clan. LD are similar to our BD. 

At the SPS energy the shoulder structure appears in MD [20"]. As it was marked in the 
branch scheme, the gluon fission is strengthened at higher energies. The independent 
evaporation of gluon sources of hadrons may be realized as single gluons as groups 
from two and more fission gluons. Following 10] we name such groups - clans. GDM 
with two kinds of clans describes data ll2lh very well (Fig. 4). Moreover "the charged 
hadron/pion" ratio at 62 GeV is equal to ~ 1.6 and agrees with Au-Au peripheral 
interactions at 200 GeV and with pp interactions at 53 GeV [18]. 

The specific feature of GDM is the dominance of active gluons in MP. We expect the 
emergence of many of them in nucleus collisions at RHIC and the formation of a new 
kind of matter (OGP) at high energy. Our QGS can be a candidate for this. 

Experiments [22] have shown that the measured cross sections of SP are several times 
larger than the expected ones from QED. The phenomenological glob model explains the 
SP excess lIlSO . We consider that at a certain moment the QGS or excited new hadrons 
may set in an almost equilibrium state during a short period of time. That is why, to 
describe massless photons, we have used the black body emission spectrum |23]. At 70 
GeV/c an inelastic cross section is equal to ~ AOmb, the SP formation cross section is 
about 4-mb \ and since Oy ~ ny(r) • a,„ then ny ^ 0.1. If Uy and temperature T{p) 
(p-momentum) are known, then we can estimate the emission region size L of SP. The 
obtained values L ~ 4 — 6/m [5, 24]. 

In conclusion one can show how GDM may explain experimental pp annihilation data 
at tens GeV/c lEsll . The differences between pp and pp inelastic topological cross sec- 
tions {Ao„{pp — pp) = On{pp) — o„{pp)) point out the contribution of different mech- 
anisms in MP. The negative values of /2 indicate the predominance of the hadronization 
stage in MP. According to MGD, the active gluons are a basic source of secondary 
hadrons. 

There are three valent ^^-pairs at the initial stage of annihilation . They can turn to 
the "leading" mesons which consist of: a) valent quarks or b) valent and sea quarks 
1 2411 . In the case a) three "leading" neutral pions (the "0"-topology) or two charged and 
one neutral "leading" mesons ("2"-topology) may form. In b) case the "4"- and "6"- 
topologies are realized. The neutron and antineutron formation is possible, too. 

The active gluons emerge together with the formation of intermediate topology. At 
this region hadronization dominates since is negative. In the simple case for m active 
gluons GF Q^{z) = [1 +n''/N{z-l)]"'^, and /2 = -m{n^'f/N. If m grows while 
increasing the energy of the colliding particles, then /2 will decrease almost linearly 



from m that agrees with the data ll25ll . 

According to GDM |5] and taking into account an intermediate charged topology 
("0", "2" and "4") with active gluons, GF Q{z) for final annihilation MD (AOnip'p - 
PP) I Ln -^cTaj {pT ~PP)) m^y be written as the convolution gluon and hadron components 

M() M2 M4 

m m m 

The parameters of cq, C2 and C4 are parts of intermediate topologies. GDM describes data 
(Fig. 6) with the ratio cq : C2 : C4 = 15 : 40 : 0.05 and Mq M2 1 - 2, M4 ~ 4 [5]. 
The carried out investigations allow one to understand deeper the MP nature. 

I thank the Organizing committee for a wonderful possibility to take part in 
ISMD2005. 
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